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SYSTEMS  OF  HYPERBOLIC  DIFFERENTIAL  EQUATIONS  Of  THE  FIRST  ORDER.  II. 


Wolfgang  Haack  and  Guenter  Hellwig 


Berlin-Cha  rlott  enb  urg 


! 

INTRODUCTION 


The  following  studies  are  a direct  continuation  of  our  first 
report  on  this  topic  which  appeared  recently  it;  this  journal.  In 
Section  I of  our  first  report  we  proceeded  ircir  a hyperbolic  system 
cf  linear  differential  equations; 

I 

I 

II* Ut  />'r.  f r«  . + r ^ 0. 

11* + V cj  + r 1/  f . ,,  I) 

(suiiimnied  over  k = 1,  2) 
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1 
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Through  invariant  derivatives  in  the  characteristic  directions 

(2) 

and  appropriate  linear  combinations 

(3) 

I rm  ti  * if  r\  » »».  M f rJ  f 

the  system  (1)  was  reduced  to  a normal  forir.  1 1;  sections  II  and  III 
existence  theorems  were  proved  for  the  functions  U,  V and  u,  v- 

Below  we  have  continued  the  numbering  of  the  sections  and 
formulas  and  in  section  IV  we  shall  deal  with  several  integral 
theorems.  Two  Pfaffian  forms  u,,  (I.li)  telong  to  the  vector 

system  (2),  Now  if  II,  V,  respectively,  u,  v are  solutions  of  (1)  and 
W,  Z are  arbitrary  functions,  then  through  the  extrinsic  derivative 
of  the  Pfaffian  form  ii  iM.  + zrn,  one  arrives  at  integral  relations  for 
the  functions  II,  V.  The  requirement  that  cortair.  terms  in  the 
integrals  vanish  loads  to  the  adjoint  system  of  differential 
equations  for  W and  Z. 


In  a manner  similar  to  that  in  the  Riemann  method  which  is  known 
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from  the  theory  of  a partial  differential  equation  of  the  second 
order,  one  can  solve  the  Cauchy  initial  value  problen  of  system  (1) 
usinq  quadratures  if  the  characteristic  initial  value  problem  is 
solved  for  the  adjoint  system.  The  result  is  the  integral 
representation  of  the  sought  functions  (IV-  28,  31).  Section  IV 
concludes  with  the  proof  that  the  classical  Kieirann  method  is 
contained  in  these  integral  representations  as  a special  case. 


In  section  V the  considerations  will  Le  transformed  for 
quasilinear  systems.  Systems  also  result  in  the  character  is»- ic  form 
which  was  investigated  by  P-  Courant  et  al.  (sec  the  first  report, 
fcotnote  [ ^ ]) . We  show  that  here  also  a normal  form  is  always  present 
and  indicate  when,  through  transposition  of  the  dependent  and 
independent  variables  a "linearization"  is  possible  and  outline 
briefly  a difference  method  for  the  approximate  solution  of  various 

initial  value  and  boundary  value  problems.  In  conclusion  as  an 
example  we  shall  deal  with  t)ie  the  plane  and  dynamically  balanced 
stationary  flow  of  an  ideal  compressible  fluid.  We  calculate  the" 
normal  form  of  the  problem  and  arrive  at  a methed  of  characteristics 
for  approximate  determination  of  the  solution  which  also  uses  the 
epicycloid  networlrs  and  which  therefore  may  offjtr  certain 
simplifications  as  opposed  to  the  familiar  methed  (cf.  sec^^ion  V, 

footnote  [ ’ T)  . 
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Finally  we  la  not  want  to  neglect  thanking  Mr.  L-  Bieberbach  for 
revision  or  the  tnanuscripts  of  both  parts  and  fcr  numerous 
suggestions  for  improvement. 

SECTION  IV. 

QUASI-niEM ANN  METHOD 

1.  Integral  theorem  of  the  normal  form.  The  construction  of  the 
characteristic  theory  of  a system  (1)  using  Pfaffian  forms,  as  was 
done  in  Section  I brings  the  advantage  that  one  can  make  direct 
statements  about  the  behavior  of  the  sought  functions  in  general 
terms.  Along  with  E.  Cartan  we  designate  the  "extrinsic  product"  of 
two  Pfaffian  forms  wi,  up  with  [ui,  u?  ] and  the  "extrinsic 
derivative"  of  a Pfaffian  form  with  [ d u ].  It  0 is  a range  in  which 
the  coefficients  of  the  form  u are  continuously  differentiable  tnen 
the  integral  theorem  is  valid  • 


fFCOTNOTE:  Cf.,  for  example,  B.  Blaschke,  "Introduction  to 
differential  geometry"  Berlin  1950.  ENC  FOCTNOTF] 


(IV.  1) 
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According  to  Section  I two  Pfaffian  forms  wi,  uj  belong  to  a 
system  of  hyperbolic  differential  oguations  (1).  They  form  the  basis 
of  a form  ring,  which  using  two  auxiliary  functions  W(x*,  x^)  , Z(xi, 
X*)  we  write  in  the  form 

( i 2 ) 4.1  =<!  II  r 1.1,  z / 1 1 

In  this  case  U,  V are  to  be  considered  as  solutions  of  the  normal 
form  (1.9) 

V.  — ^V  + BV  + C. 


The  application  of  the  integral  theorem  {IV,  1)  to  the  form  u (IV. 2) 
yields 


(IV.  4)  + J J\a<B’Vm,  + ZVm,)\. 


For  transformation  of  the  right  sid<  we  use  the  fundamental  formulas 
known  for  extrinsic  derivatives  [ * ]; 
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and 


u. 


If  wt‘  proceoii  with  this  in  (IV. M)  t lu'ii  the  inttqial  tfieoreni  becomes 


(IV. 10) 


^ I II  r M,  + z i'u,\ 


»r  rs-  zch;  iu.u. 


If  we  now  note  that  U,  V are  considered  -is  solutions  of  the  normal 
form  (IV.  1)  then  we  can  replace  and  in  (IV.  10)  according  to 

(I  V.  1)  and  o ht  a in 

/(HT..  4 ZVm,)  ^ j j r 111;+  H’(«  + S,-  ZH* !«.«.! 

Z.-Z|.<  • Hl+W.il'K-,1 
+ 11  [H  r - zc;  iw  ...!, 

This  is  a general  integral  thooitm  whicn  must  le  satisfied  by  evei y 
solution  system  of  (IV. 1)  with  any  tunctions  w , Z. 

2.  dene ra 1 iza t ion  of  tne  Kiemann  method.  It  is  easy,  usinq  the 
inteqral  theorem  (IV.  11)  to  reduce  the  taucliy  initial  value  problem 

of  (IV. d)  for  a curve  K a "char  acter  ist  ic  initial  value  problem" 
(Section  III). 


For  the  arbitrary  functions  w,  2 we  make  the  deirand 

11, - + II  i«  t o 

(IV.  12)  z.  a. 

That  is  a system  of  partial  differential  equations  of  the  hyperbolic 
type  for  the  functions  W and  Z in  the  rcrmal  form  with  the  same 
characteristics  «,  i.  We  call  it  the  adjoint  system  to  (IV.  3).  The 

adjoint  system  is  simpler  than  the  original  (IV. 3)  in  that  it  is 
always  homogeneous:  it  is  solvable,  however,  only  through  quadratures 

according  to  the  definition  of  Section  1 when  the  initial  system 
(IV,  3)  is  solvable  through  quadratures.  If  in  (IV.  1 1)  W and  Z are 
known  adjoint  functions,  i.e,,  solutions  of  (IV. 12)  then  only  known 
functions  appear  under  the  remaining  double  integral.  It  will  now 
come  to  transforming  the  boundary  integral  through  the  choice  of  an 
appropriate  range  so  that  the  sought  functions  U,  V appear  before  the 
integral.  In  the  case  of  the  Cauchy  initial  value  problem  U,  V are 
given  along  a curve  h,  the  function  values  U(P),  (VP)  are  sought  in 
any  point  P outside  of  K. 

He  consider  a region  G whose  boundary  is  termed  by  the  initial 
curve  K and  the  characteristics  through  point  F (Fig.  1).  For  this 
region  G the  prerequisites  of  the  integral  theorem  arc  satisfied.  As 
can  be  seen  from  Fig.  1 one  can  split  the  boundary  integral  up  into 
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f 
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IIV.  13) 


The  first  integral  is  to  be  taken  along  the  initial  curve  K, 
U,  V ate  known.  For  transtormat ion  of  the  integral  over  the 
characteristics  wa  consider  the  equation  (IV. 3)  and  write  it 
form 


on  which 


in  the 


(IV.  1i*) 


flF  - I’. -All  -c;  Ai'  — y--  hy-i'. 


In  order  not  to  have  to  exclude  the  zero  places  of  the  given 
functions  B and  a in  the  future,  we  want  to  write  the  adjoint/*,/* 
function  W on  the  boundary  PiP  and  the  ad  ;;cint  function  2 on  the 
boundary  PjP  in  the  fora  ^ 

,|V  15)  ” »"f  '■/*.:  z --  d / "id  /*/*, 

[FOOTNOTE;  ^ The  next  thing  would  be  to  replace  the  function  U 
according  to  (IV.  14)  by  » 1 1 v;- s r-F)  ■ in  the  characteristic  parts  of 

the  boundary  integral  (IV.  14),  respectively,  (IV.  13)  for  example  in 
the  integral  j and  then  through  partial  integration  with  respect 

to  V;  to  move  the  function  V in  front  of  the  integral.  This  can  only 
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fc€  done,  however,  if  . Thiouyh  i iitr  cduct  icn  of  the  auxiliary 

functions  • and  * according  to  (IV. 15)  this  difficulty  is 

eliminated.  END  FOOTNOTE]  i 

The  possibility  of  such  a formulation  results  directly  from  the  fact 
that  for  the  ad^joint  functions  *i , Z ch  dract  tr  ist  ic  initial  values  ate 
to  be  sought,  i-e. , we  can  have  at  our  drspcsal  W and  Z on  the 
corresponding  characteristics  through  F.  with  (IV- 1U,  15)  we  go  into 

the  subintegrals  of  (IV. 13)  and  observe  that  up  = 0 along  the 
characteristic  PiP  and  uj  = 0 along  P2F-.  Then  from  (IV.  11)  it  follows 
that 


(IV.  16) 


JtH  I m,  + /t  u,,  f ' — ( 'till 

+ fsstu'.  - .1  f - ( |„  = f flUC  zcilio.wl 

p>-,  •'  •’ 


As  a result  of  partial  integration  the  following  formulas  result: 


(IV.  17)  ^ di' »■,  j I'ld.  ♦ wd"-*.- 

' ^ ^ /*,  < I*  't 


|iv.  18)  j»ii.  41  /r,®,  + /»(■«), 

'•6  It-,  It, 


J 


t)(H'  i'AUK  11 

Now  tor  •,,!<  <1  loll. I till'  roi  1 1 spiMi.l  i n-j  i-li  .it  ii  i t i>i  i :;f  if;:  thioiiiih  P w»' 

m.ik**  the  ilenAtKl: 


( 1 V . 1 '> ) 


H:  I N't  » <•  Hilt 


U V.  .’(1) 


tt.  1 .tW  - tl  Kllf 


Ihorrfrom  tor  .S  in  .i  i .iiiiUnii  point  P**  ol  11\,  .iii.l  loi  W in  ,i 
innilom  point  P*  of  PP,  (;.oo  Kii].  1) 


(IV,  ^'1) 


r»* 

/ • 


(I 


.ml  rr, 


I <■ 

Wi/“i-  .ml  rr, 


In  t h i r.  o.»  so  t' , .iiul  r.  .no  i nt  «'.|i  .it  i on  toiisl.iiits  whioh  wo  sh.ill  I.ifoi 
ha  VO  .1 1 on  t .1  i ;..i'oi.  .t  1 . 


'thon  in  .»<■<•<)  i d.i  in'o  with  (IV.1‘>),  foi  W,  Z lillows  t lio  iloraan.i: 


(I  V..’  I) 


/-«,  ( (./~ 


/ * 

I . 


.ml  rr 


L 


(I  V-  ^’4) 


II  r I II  I- 


.11'  rr 
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These  are  characteristic  initial  values  accordinq  to  Section  III  for 
the  solutions  W,  Z o£  the  adjoint  system,  Accordinq  to  these 
stipulations  the  integral  theorem  (IV.  16)  finally  acquires  the  form: 


« /'.r  /-I 


(IV.  25) 


,(  /■  I i/'.  = r ( /;,  _ y.f.  . 

- Ii^(  .1,  . /.;(  ,,,  . / II  I * X / « > 

' I-,  r,  I , 

- J I ;ii  < -/(  ; u,  //ill  X' 


We  want  to  summarize  the  result  up  to  this  point  in  the  following 
manner: 

If  W,  Z are  solutions  of  the  adjoint  system  (IV.  12)  which 
satisfy  the  characteristic  initial  conditicns  (IV. 2J),  (IV. 24)  with 
arbitrary  constants  Cj,  c^  then  the  linear  combination  .51  in  the 

point  P can  be  expressed  through  the  initial  values  and  known  values 
in  accordance  with  (IV. 25).  Here  #<,.•<  are  explained  by  (IV. 21), 

(IV. 22).  We  want  to  call  w,  z g uasi-Riemann  functions  and  their 
existence  is  assured  according  to  Section  11  and  III. 

A further  goal  will  be  to  obtain,  through  appropriate  choice  of 
the  constants  Cj,  c^,  two  independent  linear  combinations  from  which 
0,  V can  be  calculated  in  point  P.  Tliis  can  be  done  in  various  ways. 
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If  one  ilosires  certain  symmetry  in  the  foiirulas,  then  through 
the  demand 

(IV.  26)  I.  1;  .>)(/•)«  1 

one  can  uniguoly  determine  the  constants  C,  and  in  (IV.  21,  22)  and 

the  functions  K,  Z according  to  (IV-2J,  24).  A second  function  pair 

tt,  ^ we  subject  to  the  demand 

(IV. 27)  11.  w./'i-i; 

and  correspondingly  obtain  uniguely  determined  functions  u . . If  we 
regard  the  abbreviation  H{W,  Z)  given  in  {IV. ZS)  then  it  folLous 
directly  that 

(I  V.  2«a)  ' ~ ! !«'»  . zi  . >h».  Zi). 

(IV.  28b)  i i/’i  ,i:;/iiT,zi  ihu  . /y,. 

If  the  characteristic  initial  value  problcB  is  solved  tor  the 
adjoint  normal  form  then  the  solution  of  the  Cauchy  initial  value 
problem  of  the  system  (IV.  J)  toi  every  initial  curve  K can  be 
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cepresent^ed  by  quad  rat  ulos. 

Shorter  integral  representations  can  Le  ottained  in  the 
following  manner:  We  choose  the  integration  constants  in  (IV.  21,  22) 

so  that 


(IV-  29) 


1 Si*  l’>  I : d II  .iiif  /■/',. 


(IV. 30) 


" I Sll  „ a,.f  /■/■, 


Thereby  W,  Z,  W,  Z are  uniquely  determined.  Froir  (IV.  2'))  for  the 
sought  functions  IMP),  V (P)  it  follows  immediately  that: 


(IV.  31a)  f ,/*.  = ly.i-i  jSiW  i.,  . /■  ill  I u,  + z(  «i,i 

I'.C 

I I ( III  r \ 


(IV.  31b)  I Ti  ,V/',i  I i/',i+  id'll*  /'  ill  Tk,  t Z'lii.i 

♦ j I in  r Z(  I ,w,ia  |. 

From  the  integral  representation  (IV. 31)  it  follows  particularly 
graphically  that  tl,  V are  continuously  dependent  on  the  initial 
values.  Finally  let  us  note  that  integral  representations  for  the 
functions  u,  v of  the  Cauchy  initial  value  problem  of  the  general 
system  (1)  are  given  directly  by  the  integral  theorems  (IV.2rt), 
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(IV. 31). 


3.  Relations  to  the  "Piemann  method."  If  a -general  hyperbolic 
differential  equation  of  the  second  order  is  given  in  the  form 

(1.17),  in  which  we  designate  the  independent  variables  with  x,  y 


(IV.  32) 


t •/.  + 6«,  * r#  I / -•  0, 


then  through  the  introduction  of  both  functions 


(IV. 33)  .=1. 


we  can  write  these  as  a system.  According  to  (1,19)  in  addition  there 
results  the  normal  form 


(IV.  34) 


r.  - -fcr+  r. 

r,  - ih,  + ah  - r)l' - a V - f 


with 


>1  = 1’;  h»  + r (•; 


since  according  to  (1.18)  and  (1.1)  r.  i,. 


Two  Pfaffian  forms  (1.13)  belong  to  (IV. 34) 
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M,  •“  nV;  ■«,  i/y 

Frcin  the  integral  theorem  (IV-25)  we  immeil  iatel  y obtain  infinitely 
many  integral  representations  for  a Cauchy  initial  value  problem  from 
(IV. 32).  It  is  now  cf  interest  whether,  and  in  the  given  case, 
through  which  ''mark.ed''  demands  these  integral  representations 
directly  yield  the  familiar  "Piemann  scluticn  formula. " 

I 

He  shall  show:  If  for  the  characteristic  initial  value  problem 
cf  the  "adjoint  system"  (IV. 12),  in  addition  to  requirements  (IV. 19, 
20)  we  make  the  demand  (IV. 29)  for  the  determination  of  the 
integration  constants  then  (IV. 31a)  yields  precisely  the  Piemann 
solution  formula. 

The  adjoint  system  (IV.  12)  becomes 

(I  V.  36)  II . II  -I  /.  II. 

Z,  Z/l  - II  t •!><  - I ' II 

The  requirements  for  the  characteristic  initial  value  problem  cf 
(IV.  36)  according  to  (IV.  15)  and  (IV. 20)  since  E = 1 (IV.  37)  ii  - bW 
= 0 on  PP,  with  W(P)  = 1 

(IV.  37a)  Z = 0 on  PP^  according  to  (IV. 29). 


r 
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From  (IV.  36)  wo  eliminate  7,  by  differentiating  the  first 
equation  with  respect  to  x and  introduce  /..  ftcm  the  second 
eq  uat ion. 

Following  simple  transformation  it  follows  that 

(I  V.  38  ) II  - 'll  11;  - II  /'  , • II  < ” 

(IV. 38)  is  precisely  the  "adjoint  hiemann  equation"  which 
belongs  to  (IV. 32) . 

According  to  (IV.  36)  the  requirement  (IV-  17a)  goes  over  into 

(IV.  39)  H,  - aW  = 0 on  PP^. 

(IV. 37,  39)  are  the  "Piemann  demands"  which  belong  to  (IV. 38) 

According  to  (IV. 33,  34)  we  set 

(IV. 40)  ' - f . I 


Now  we  arrive  directly  at  the  Rioirann  solution  formula  if  we 


DOC  = 0926 


PACE  18 


substitute  the  values  (IV. 40)  into  (IV.JIa).  Considering  the  first 
equat  ion  in  (I  V.  J6  ) 

tin  - M /■{»’(*»  + fa) »'«)»rf|r> 

(IV.41)  f.V 

- / ffWHiilyl 


This  is  the  Hiemann  representation  formula  for  the  Cauchy  problem  of 
the  hyperbolic  differential  equation  (IV- J^:).  (Cf.  perhaps  Goursat 
"Cours  D'Analyse"  Paris  1927,  Volume  3,  page  150- ) 

SECTION  V.  SYSTEMS  OF  gUASILINEAH  DIFFEBENTIAL  EQUATIONS. 

In  the  introduction  it  was  already  pointed  out  that  in  the  past 
ten  years  in  connection  with  certain  prcblems  cf  fluid  dynamics 
methods  of  characteristics  were  developed  for  the  solution  of 
hyperbolic  systems  of  quasilinear  differential  equations  with  two 
independent  variables 


[FCOTNOTE;  ■>  Even  today  t)ie  abundant  literature  is  accesible  here 
only  with  difficuly  so  that  we  can  only  mention  the  follow inq: 

P.  Sauer,  "Theoretical  Introduction  in  Gas  Cynairics,"  Berlin  1 94  3. 
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Kl,  Oswatitsch,  "flt>thodii  ot  characteristics  of  Hydromechanics,"  ZAMM 
1947,  Issues  7,  8,  9. 

Courant-Fr iedrichs , "Supersonic  Flow  and  Shock  Waves,"  New  York  1948. 

A.  Ferri,  "Elements  of  Aerodynamics  of  Supersonic  Flows,"  New  York 

1949.  I 

! 

Abundant  b iblioyra ph ica 1 information  can  be  found  in  the  above  works. 

END  FCOTNOTE]  ' 

In  particular  R.  Courant,  K.  0.  Friedrichs  and  Lax  have  recently 
studied  the  quasil inear  systems  in  detail  and  given  the  prerequisites 
fo<*the  existence  of  unique  solutions  of  the  Cauchy  initial  value 
problem  (cf.  Report  I.,  footnote  2).  Here  we  want  to  limit  ourselves 
to  pointing  out  the  connection  of  the  quasilinear  systems  with  the 
theory  of  the  linear  systems  developed  in  Sections  I through  IV  and 
to  illustrating  the  advantages  using  a simple  example  of  the  plane 
and  dynamically  balanced  flow  of  compressible  liquids. 

1«  Quasilinear  systems  in  the  "characteristic  form." 

let 
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(V-1)  ■*«.  . ^i-.  + r — 0.  (summed  over  ic  = 1,  2). 

'»•  M.  i r —0 

be  a quasilinear  system  in  wiiich  «*.»<•...<•  are  functions  of  the 
independent  variables  x>,  and  of  the  sought  solutions  u,  v.  If  we 
follow  the  developments  of  Section  I and  assume  that  a solution  u,  v 
is  known,  then  the  system  (V. 1)  is  called  hyperbolic  with  respect  to 
this  solution  if  the  determinant  (1-3)  has  exactly  two  real  roots 

Therooy  the  characteristic  directions  became  dependent 
on  the  solutions  u,  v,  however,  equations  (1.  1-6)  of  Section  I 
maintain  their  full  validity,  while  the  factors  pj,  pj,  *4, 
now  be  functions  of  x‘,  x^,  u,  v.  In  this  manner  we  always  arrive  at 
the  following  "characteristic  system”  in  which  we  use  the  indices  a, 
a instead  of  points  in  the  differential  equations  (1.11)  of  the 
character! sties: 

III  « i'. 


(V.2a) 


II  y. 

y. 


The  are  known  functions  of  x*,  x^,  u,  v.  In  addition,  according 

to  (16)  there  are  the  accompanying  differential  equations  of  the 
functions  u,  v 


The  differential  equations  (V, 2d)  of  the  characteristic  curves  can 
also  be  written  in  the  form: 

(V.2c)  ll  a' — iiV;  = il  II  1 ' r’  It 

Now  the  equations  are  hocoqeneous  and  are  f;ies€Lved  during  the 
introduction  of  integrating  factors-  Hence  the  directional 
derivatives  can  be  interpreted  as  partial  derivatives.  (V. 2b)  and 
(V,2c)  form  a system  of  four  partial  differential  equations  for  the 

functions  x*,  x^,  u,  v of  two  characteristic  variables-  The 
properties  and  the  questions  of  existence  for  such  systems  are 
studied  by  R.  Courant  et  al.  (cf.  fcctrcte  I^). 

2.  Normal  form  of  quasilinoar  systems- 

Our  theory  of  linear  systems  is  based  essentially  on  the  normal 
form  (see  Section  I).  We  turn  to  the  question  of  whether  the  system 
of  equations  (V.2b)  possesses  a normal  form  in  the  sense  of  (19).  One 
immediately  recognizes  the  correctness  cf  the  theorem: 
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a)  If  in  (V.  2b)  are  functions  of  x‘,  alone,  then  one 

arrives  at  the  noraal  form  ty  sottinq: 


{V.  J) 


t'  “ 9,  " ^ I «-  «,  « t «.  > 


The  normal  form  rea^ls 


(V.4) 


r.  (;  f;  - ( 


In  this  case  C and  C are  known  functions  of  x>,  x^,  U,  V. 


Mr.  J_  Nitsche  pointed  out  to  us  that  a correspond inu  theorem  is 


also  valid  in  the  general  case' 


b)  If  v«.«4  in  (V.  2b)  arc  differentiable  functions  of  u,  v,  *•, 
then  the  system  possesses  a noiiral  form. 


We  have  to  show  that  instead  of  u,  v,  new  functions  U,  V of  u, 

V,  r can  be  introduced  which  reduce  (V.  21)  tc  tne  form  (V-4)  . For 
this  purpose  we  multiply  the  first  equation  (V.2b)  with  a function 
M(u,  T,  j‘l  and  the  second  equation  with  N’  (u,  v,  y)  and  define  M and 
N so  that  the  fallowing  equations  become  possible: 
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(V.5) 


Vti.M.  + — V,  ♦ itn.  r,  x'), 

■V  + .V  «,  fi  i=.  I-';  + i («,  r,  X*). 


That,  is  cettdiiily  the  case  if  K and  N satisfy  the  equations: 


(V-6) 


< -o.  0. 


In  these  equations  the  r are  to  be  interpreted  as  constant 
parameters.  Then 

('  •«  ♦iM.  r.  x*!;  r Vim.  r.  x‘) 


such  that 


(V.6a) 


r.  - - If,., 

r.  V. -.v«,. 


= If,, 

r,  — V.  -- .%  «, 


Ey  differentiation  in  the  direction  of  the  characteristics: 

f.  ■=•  + .V,,  i 0„a‘ 

y-  =.  .\  «,  M;  V «,  1;  + V,  i' . 


Here  the  last  terms  are  known  functions  of  u,  v,  r.  if  the 
integrating  factors  M,  N are  defined  according  to  (V.6).  The 


func^’iondl  determinant  of  the  functions  U,  V with  respect  to  u,  v is 


NN  (p,#2-p2ff,)  . Cn  account  of  the  independence  cf  ».  • , d,  V are  also 

independent  with  respect  to  u,  v-  Hence  it  follows  inaiediately  that 
(V.2b)  can  be  written  in  the  normal  form  (V-4). 

c)  Kor  a quasilinear  system  in  the  noittal  form  (V.2a,  4)  the 

integral  theorem  ( IV.  1 1)  can  be  accepted,  lor  any  two  solutions  U,  V 
of  the  guasilinear  system  and  two  arbitrary  functions  w,  Z according 
to  (IV. 11)  the  integral  theorem  is  valid: 


y'lir  I u,  zr*.,  =.  I fy  \ it;  + » x | iu.«,i 
(V.  7)  + J fi'  [ - /.  XH  \ u u, 

^ f I - ZCi  Iw.i. ,) 


whereby  R,  S are  given  by  (IV.  9). 

By  setting  the  bracl^eted  expression  in  both  double  integrals 
equal  to  zero  the  adjoint  equations  result.  If  cnc  assumes  U,  V to  be 
known  then  the  Pfaffian  forms  can  be  made  integrable  and  the  adjoint 
system  can  be  reduced  to  the  simple  form 


DDC 


0 'I (1 


I’ AUK 


air  thr  rij  uat  i on:.;  t oi  lioth  systrm;.  of  ( In  i.u:t  oi  i;?t  if;;. 

il)  Li  nrar  i zah  1 r y r.t  rms.  C t j oi  app  J ir  at  ions  is 

thr  sprcial  c.*!;*-  t oi  wliicli  in  (V.U) 

C ^ C ~ 0 

Ono  can  1 1 ansp<.);;r  tin*  ii)lr  cl  thr  iirp*‘inlrnt  and  j luioprniion  t 
variables,  i.o.,  one  i n t t>i  pit't  s li,  V .is  .in  inilrpiMuh'iit  variable  in  a 
U,  V-plane  and  considers  j (U,  V)  .is  sou.)  lit  tnnit  ions  ot  t In' 
independent  vaii.ibl.'s  U,  V.  Th«'n  (V.2c)  .if'jiear  a:;  a sy.stmi  ol 
differential  eijiMt  ions  in  c ha  i act  e r i ::t  ic  tcim  tor  the  Innctions  x'  (U, 
V)  . The  cui  ves  II  - con:;t,  V ~ I'onst  ,u  r t h*  char,u:t,ei  ist  ics  ot  the  11, 
V-plane,  which  are  lixisi,  i.e,,  .u  »>  indeptndent  ot  the  .'^oitiiht 


functions  i 

' (11 » V ) . Tlie  r r w i t li 

o no 

twi.*; 

the  s.iine  r i*  1 .i  t i ons  h i ( 

d ;i 

\ n 

the  1 inear 

h ypi'i  bo  lie  s y;.  t e ITS 

<1  iiij 

ca  11 

c.irry  the  tlu'ory  ove  i 

1 i i ect  1 y 

We  shall  be 

:;a  t i sf  i ed  with  t h i 

1 ' i.» 

h 0 L' t 

irteiiiicr  (si'e  undi’i 

(he) 

bu  t 

would  still  like  to  iiot.>  that  .ilso  in  qrnri.il  cases  such  a 
1 inear  i zat  ion  is  possible,  for  example,  if  C,  and  ' « are 

tuiK'tions  .ilone  ot  1),  j'  and  V,  ' , i I's  pcct  i vr  1 y . 


1.  Urneral  difteieiui'  iretho.l. 
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From  the  equation  system  (V.2a)  and  (V.4)  (resp.  (V,2b>)  crie 
arrives  directly  at  difference  methods  which  a p pt ox iaa t el y solve  the 
Cauchy  problem  and  all  of  the  mixed  initial  value-  anil  boundary  value 
problems  discussed  in  Section  III.  Let  us  first  assume  that  the 
solutions  of  the  systems  (V.2a)  and  (V.4)  are  known  and  s, 
respect i ve 1 y , s are  the  curve  lengths  along  the  characteci st ics.  Then 
the  equations  can  be  written  in  integral  form.  Approximately 


1.  r 

a J'i 

./ 

It  / 

/ * 

u ! >. 

n 

. J 

J.f 

i J- 

- t 

Along  the  curve  K let  the  initial  values  D (K)  , V(K)  be  given;  along 
the  boundary  R let  some  function  •i>  (u,  v),  respectively,  ((I,  V)  be 
known  (see  Fig.  2)  . Beginning  from  the  intersection  0 one  passes  the 
characteristic  directions  through  a number  of  points  (1,  2,  3)  of  K. 

They  are  known  as  quotients  of  the  first  two  equations  (V. 8)  along  K. 
The  characterist  ic  straight  lines  intersect  in  pairs  in  the  points 

6,  7.  If  one  normalizes  — I and  • l.  then  As,  respectively  As 

is  the  spacing  ot  the  points  on  the  a-,  respectively,  a-line.  Cne  can 
thus  calculate  j.v  and  JiV  and  determine  1),  V in  the  points  “j,  6, 

7. 
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Assuminvj  tliat  the  chat  act  cc  is  t ic  approaching  in  boundary  point  4 
belongs  to  system  II  of  the  equations  (V. 8)  then  in  4 one  knows  the 
function  V and  from  V and  the  boundary  condition  vlU,  V)  can 

calculate  the  function  U (cesp.  u,  v) . In  this  case  it  is  assumed 
that  * is  uniquely  solvable  with  respect  to  U.  The  method  can  he 
continued  in  the  same  manner  in  the  range  lying  between  K and  R.  Uith 
respect  to  the  boundary  curve  R it  is  to  be  assumed  that  in  the  range 
between  K and  P there  is  one  and  only  one  characteristic  passing 
through  0, 

4)  steady  flow. 


As  an  example  dynamically  balanced  and  plane  steady  flow  c£  an 
ideal  gas  will  be  treatCvl  in  scmewhat  mere  detail.  In  view  of  the 
abundant  literature  on  the  subject  we  may  be  brief.  In  cylinder 

coordinates  r,  x let  #(r,  x)  be  the  velocity  potential  and  a be  the 
local  velocity  of  sound.  Then  the  following  equation  is  valid: 


(V-9) 


v..  1 1 


1 I, 


V 


In  this  case  with  the  constants  * i.  «•.  » 


a* 


(V-10) 


• 1 


(wj—  #r*i,  «■’  gj. 


then  from  (V.9)  and  the  condition 

^ ^ r 

(V.  12)  I «•••<•»  »! 


there  results  a quasilinedr  system  of  differential  equations  foi.  the 
functions  w,  9,  which  qive  the  contribution  and  the  direction  of  the 
flow  velocity  toward  the  x~axis.  For  w > a the  systea  is  hyperbolic. 


I 


F 


Here  it  is  recoaaended  to  not  go  directly  to  the 


"characteristic"  derivatives,  a,  a but  to  gc  first  to  another  systen 
cf  "invariant"  derivatives. 


There  are 


i M .r  r<M  # 
r tr  fiiii  . 


(V.13) 


— «■  Ip  y dill  ♦ 
tr  (jT  j*  con  ^ 


two  vectors  of  whicn  (3  coincides  with  the  velocity  ">  and  (3  is 
perpendicular  to  it.  In  this  case  y is  the  so-called  "Hach  angle" 


11 
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with  sin  j = a/w  (cf.  Fiq-  i)  . The  invaiiant  derivatives  of  a 
function  S in  these  directions  are 


If  one  calculates  these  derivatives  for  the  functions  w and  8 and 
observes  (V,4)  and  (V.  12)  then  the  system  of  invariant  differential 
equations  follows  after  a simple  intermediate  calculation; 


(V.  15) 


„ n-;-  i., 


ft 


II 


a)  Natural  equations.  The  formulas  (V15)  allow  a simple 
geometric  interpretation.  It  one  designates  the  curve  length  with  s 
and  the  bend  radius  of  the  flow  lines  with  <’■  then  on  account  of 

(V. 14) 


(V.  16a) 


•ID  I 


I' 


/ • I 

IW.-yni  I 

•**  f. 


If  WG  call  n and  a,  the  curve  length  ami  the  bend  radius  of  the 


crthogonal  trajectories  of  the  flow  lines,  then  according  to  (V.l4) 


i 
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Ey  introduction  in  (V.  15)  the  so-called  natural  equations  ot  the  flow 

.<  In  » . , ,1  HU  • 

J.  - r I 


{V.  17) 


.ill.  » 
Urn 


F‘'3  3 

They  can  be  obtained  quite  effortlessly  in  this  manner-  One 
recognizes  then  that  in  the  case  of  tlcw  along  a wall  w only  changes 
continuously  if  the  curvature  of  the  wall  is  ccntinuous. 


b)  The  normal  form  ot  the  system.  In  order  to  reduce  the  system 
of  equations  to  the  characteristic  form  we  have  to  add  both  equations 
(V. 15)  once  and  subtract  them  once  and  finally  divide  them  by  w.  If 
in  addition  we  set 


then  we  ob*^ain  the  characteristic  system 


(V.18) 
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Here  the  c hdract.er  i st  ic  diretrtions  a,  a,  as  a comparison  with  (V-14) 
shows,  are  qiven  by 


(V.  19)  ■■  ' " • 

. ■'  Hi  it 

Since  in  (V.  18)  the  factor  is  a function  of  w alone  the  system 

possesses  a normal  form  with  the  functions 

(V.20)  ' / 'I 

The  normal  form  reads: 

(V.21)  ' I ■ I 

where  w,  9 are  to  be  considered  as  functions  of  U,  V. 

e)  Let  us  next  examine  the  limitinq  case  ■< 

Equations  (V.  19)  to  (V.21)  qo  over  into  the  formulas  tor  plane 


1. 

t 


flow  for  which: 


I 
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The  curves  U = const,  V = const  form  a fixed  ytid  in  the  w,  »-plane 
and  arc  the  characteristics  of  the  linear  system  (V.19)  if  one  writes 
it  in  accordance  with  (V.2c)  in  - . and  ■ - . The  linearization  of 

the  differential  equation  of  plane  flow  by  transposition  of  the 
dependent  and  independent  variables  has  been  kncwn  for  a lonq  time 

(cf,  perhaps  R.  Sauer  [>]  page  111).  The  curves  U = const,  V = const 
are  epicycloids  of  the  w,  S-plane.  Epicycloid  grids  are  available  in 
various  scales  and  are  widely  used  for  approximate  determination  of 
plane  flow. 

d)  Cycloid  grids  can  also  be  use  to  advantage  tor  dynamically 
balanced  flow,  which  to  the  best  of  our  knowledge  has  not  been 
mentioned  in  the  literature  up  until  new.  We  therefore  want  to 
briefly  explain  the  method: 

Along  the  characteristics  we  choose  the  curve  lengtlis  s and  s as 
parameters.  Then  (V.19)  and  (V.21)  can  te  written  approximately  in 
differences: 


/ * I ^ * I . I .»  / . 

fV.22)  ,,  J • -in  . It  ; . / t , 


If  in  two  points  P,  Q of  the  x,  r-plane  one  knows  the  values  of  w,  9 


(Fig.  4)  then  f tom  the  first  two  equations  (V.22)  one  calculates  the 
characteristic  tangents,  their  point  of  intersection  R and  obtains  As 
and  From  them  AU,  AV  can  be  determined.  The  values  of  w (R) , »(R) 

can  now  be  effortlessly  read  off  of  the  cycloid  grid.  In  order  to 
show  that,  we  proceed  from  the  fact  that  indefinite  integrals  appear 
in  (V.20).  Since  vanishes  for  w = a we  can  introduce  the 

definite  integrals  in  (V.20) 


Then  for  w = a 

•'  ' I . 

The  value  AW  appears  therefore  on  the  icner  boundary  circle  of 

the  grid  as  AU  = -A9  and  correspondingly  AV  as  AV  = ♦Ad.  How  the 
following  wethod  results  (Fig.  5): 


Fig.  4.  ((KEY:  1)  X,r-plane.)) 


Fig.  5.  ((KEY:  1)  w,»-plane.)) 


DOC 


From  the  given  w(P),  8 (P)  and  m (Q)  , • (g)  by  drawing  the 
corresponding  radius  vectors  ot  the  w, 8-plane  one  immediately  finds 
the  points  po,  QO.  On  the  U-line  of  the  cycloid  grid  one  approaches 

the  boundary  circle  (w=a)  through  go,  marks  oft  the  value  AU  = -A*, 
calculated  from  the  thiru  eguation  (V.22),  in  radian  measure  and 
turns  back  on  the  new  grid  line  D ♦ AtJ.  If  one  does  tJ.e  same  thing 
with  P,  as  a section  of  the  new  grid  lines  one  obtains  [illegible] 
can  read  off  w(R),  8(R)  directly.  Since  the  transition  froa  U to  U AU 

(resp.  from  V to  V ♦ aV)  signifies  only  a stable  rotation  of  the  grid 
curves  by  the  angle  (-AU),  resp-,  AV,  the  process  can  be  further 
techan ized . 
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